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Weak lensing reconstruction through cosmic magnification I: a 
minimal variance map reconstruction 
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ABSTRACT 

We present a concept study on weak lensing map reconstruction through the cosmic magni- 
fication effect in galaxy number density distribution. We propose a minimal variance linear 
estimator to minimize both the dominant systematical and statistical errors in the map recon- 
struction. It utilizes the distinctively different flux dependences to separate the cosmic mag- 
nification signal from the overwhelming galaxy intrinsic clustering noise. It also minimizes 
the shot noise error by an optimal weighting scheme on the galaxy number density in each 
flux bin. Our method is in principle applicable to all galaxy surveys with reasonable redshift 
information. We demonstrate its applicability against the planned Square Kilometer Array 
survey, under simplified conditions. Weak lensing maps reconstructed through our method 
are complementary to that from cosmic shear and CMB and 21cm lensing. They are useful 
for cross checking over systematical errors in weak lensing reconstruction and for improving 
cosmological constraints. 
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1 INTRODUCTION 

Weak gravitational lensing has been establishe d as one of the 
most powerful probes of the dark un i verse (Refregierj 120031: 
lAlbrechtetaHl2006| ; iMunshi etal1l2008l : iHoekstra & Jainfl2008h . 
It is rich in physics and contains tremendous information on dark 
matter, dark energy and the nature of gravity at cosmological scales. 
Its modeling is relatively clean. At the multipole £ < 2000, grav- 
ity is the dominant force shaping the weak lensing power spec- 
trum while complicated gas phys ics only affects the lensing power 
spectrum at less than 1% level Jwhitdbool ; Izhan & Knoxll2004 
lling et al.ll2006l : iRudd et alj|2008h . This makes the weak lensing 
precision modeling feasible, through high precision simulations. 

Precision weak lensing measurement is also promising. The 
most sophisticated and successful method so far is to measure the 
cosmic shear, lensing induced ga laxy shape distortion. After the 
first detections in the year 200 dBaconetal.ll2000t iKaiser et al.l 
l2000l : IVan Waerbeke et alj2000tlWittman et alfcOOOl) , data quality 
has been improved dramatically (e.g. JFu et al.|2008l) . Ongoing and 

S'aned surveys, such as DEsQ, LSST0, JDEMq and Pan-STARRS 
have great promise for further significant improvement. 
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However, weak lensing reconstruction through cosmic shear 
still suffers from practical difficulties associated with galaxy shape. 
Thes e include shape measurement errors (addit ive and multiplica- 
tive) jHevmans et al.ll2 006: Massey et al. 2 007 ) and the galaxy in- 
trinsic alignment dCroft & Metzlerll2000tlHeavens et alfcOOOl: ling 
| 2002| : |Hirata & Seliakll2004l:lMandelbaum et alJl200fj| ; lHirata et"aL 
l2007l : lokumura. & Iindl2009l : IOkumura et al.l2009l) . 

An alternative method for weak lensing reconstruction is 
through cosmic magnification, the lensing induced fluctuation in 
galaxy (or qua sar and any other celestial objects) number density 
dMenaral2002l and references therein). Since it does not involve 
galaxy shape, it automatically avoids all problems associated with 
galaxy shape. 

However, the amplitude of cosmic magnification in galaxy 
number density fluctuation is usually one or more orders of mag- 
nitude overwhelmed by the intrinsic galaxy number fluctuation 
associated with the large scale structure of the universe. Ex 
isting cosmic magnification measurements ( Scranton et al. 1200: 
Iffildebrandt et al.ll2009l ; iMenard et al1l2010l : Ivan Waerbe5l201 
circumvent this problem by cross-correlating two galaxy (quasar) 
samples widely separated in redshift. Unfortunately, the measured 
galaxy-galaxy cross correlation is often dominated by the fore- 
ground galaxy density-background cosmic magnification correla- 
tion and is hence proportional to a unknown galaxy bias of fore- 
ground galaxies. This severely limits its cosmology application. 

The intrinsic galaxy clustering and cosmic magnification have 
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different redshift and flux de pendence, whi ch can be utilized to 
extract cosmic magnification. IZhang & Pen] (2006) demonstrated 
that, by choosing (foreground and background) galaxy samples 
sufficiently bright and sufficiently far away, the measured cross 
correlation signal can be dominated by the cosmic magnification- 
cosmic magnification correlation, which is free of the unknown 
galaxy bias. However, even for those galaxy samples, the fore- 
ground galaxy density-bac kground cosmic ma gnification correla- 
tion is still non-negligible jZhang & PerJl200q) . This again limits 
its c osmology applicatio n due to the galaxy bias problem. 

IZhang & PerJ d2005h further showed that, since the galaxy in- 
trinsic clustering and cosmic magnification have distinctive depen- 
dence on galaxy flux, cosmic magnification can be extracted by 
appropriate weighting over the observed galaxy number density in 
each flux bins. Weak lensing reconstructed from spectroscopic red- 
shift surveys such as the square kilometer array (SKA) in this way 
can achieve accuracy comparable to that of cosmic shear of stage 
IV projects. These works demonstrate the great potential of cosmic 
magnification as a tool of precision weak lensing reconstruction. 
Furthermore, since cosmic shear and cosmic magnification are in- 
dependent methods, they would provide valuable cross-check of 
systematical errors in weak lensing measurement and useful infor- 
matio n on galaxy physics such as galaxy intrinsic alignment. 

IZhang & PerJ J2005t) only discussed two limiting cases, com- 
pletely deterministic and completely stochastic biases with known 
flux dependence. In reality, galaxy bias could be partly stochastic. 
Furthermore, the flux dependence of galaxy bias is not given a pri- 
ori. In this paper, we aim to investigate a key question. Are we 
able to simultaneously measure both cosmic magnification and the 
intrinsic galaxy clustering, given the galaxy number density mea- 
surements in flux and redshift bins? 

The paper is organized as follows. In Sj2] we present the ba- 
sics of cosmic magnification and derive a minimal variance esti- 
mator for weak lensing reconstruction through cosmic magnifica- 
tion. In a companion paper we will discuss an alternative method 
to measure the lensing power spectrum through cosmic magnifi- 
cation. In Sj3] we discuss various statistical and systematical er- 
rors associated with this reconstruction. In |4] we target at SKA 
to demonstrate the performance of the proposed estimator. We dis- 
cuss and summarize in ij5] We SKA specifications are specified in 
the appendix [A] Throughout the paper, we adopt the WMAP five 
year data dKomatsu etai]|2009l) with Sl m = 0.26, Q a = 0.74, 
Q b = 0.044, h = 0.72, n s = 0.96 and ct 8 = 0.80. 



2 A MINIMAL VARIANCE LINEAR ESTIMATOR 

Weak lensing changes the number density of background objects, 
which is called cosmic magnification. It involves two competing 
effects. A magnification of solid angle of source objects, leading to 
dilution of the source number density, and a magnification of the 
flux, making objects brighter. Let n(s,z s ) be the unlensed num- 
ber density at flux s and redshift z s , and the corresponding lensed 
quantity be n L (s L , z B ). Throughout the paper the superscript "L" 
denotes the lensed quantity. The galaxy number conservation then 
relates the two by 

n L (s L , z s )ds h = —n(s, z s )ds , (1) 



where s L = sfi. n(0,z s ) is the lensing magnification at corre- 
sponding direction 6 and redshift z s , 

At(0,z s ) = - 3_ ~ 1 + 2k(0,z s ). (2) 

(l- K (0,z s )) 2 -7 2 (0,^) 

The last expression has adopted the weak lensing approximation 
such that the lensing convergence k and the lensing shear 7 are both 
much smaller than unity (jftj, I7I <C 1). ft for a source at redshift 
z s is related to the matter overdensity S m along the line of sight by 

,„ . 3H§Sl m f Xs D( Xa - x)D(x) , , , s, 

*) - J o % (Xs) M*- *)(i + *)dx ■ 

(3) 

Here, \ an d Xs are the radial comoving distances to the lens at 
redshift z and the source at redshift z s , respectively. D(x) denotes 
the comoving angular diameter distance, which is equaling to \ f° r 
a flat universe. 

Taylor expanding Eq.[T]around the flux s L , we obtain 

n L (s L ) = \n{ — ) « n(s L )[l + 2(q - l) K + 0(ft 2 )] , (4) 

Where the parameter a is related to the logarithmic slope of the 
luminosity function, 

_ s L dn(s L ) ^_ dlnn(s L ) ^ 
~ n(s L ) ds L dlns L 

Notice that ra(s L ) is related to the cumulative luminosity function 
N(> s) by N(> s) = n(s L )ds L . Weak lensing then modifies 
the galaxy number over-density to 

5g ~ 5 S + 2(a — 1)k = S g + gn . (6) 

Here, <5 g is the intrinsic (unlensed) galaxy number over-density 
(galaxy intrinsic clustering). For brevity, we have defined g = 
2(a— 1) and will use this notation throughout the paper. Obviously, 
to the first order of gravitational lensing, the cosmic magnification 
effect can be totally described by the cosmic convergence and the 
slope of galaxy number count. 

The luminosity function averaged over sufficiently large sky 
area is essentially unchanged by lensing, 

{n^ = {n)(l + 0({K 2 )))~{n). (7) 

The last expression is accurate to 0.1% since (ft 2 ) ~ 10 -3 . The 
above approximation is important in cosmic magnification. It al- 
lows us to calculate a by replacing the unlensed (and hence un- 
known) luminosity function n with n L , the directly observed one. 
This means that a is also an observable and its flux dependence is 
directly given by observations. We will see that this is the key to 
extract cosmic magnification from galaxy number density distribu- 
tion. 

The biggest challenge in weak lensing reconstruction through 
cosmic magnification is to remove the galaxy intrinsic clustering 
<5g. This kind of noise in cosmic magnification measurement is 
analogous to the galaxy intrinsic alignment in cosmic shear mea- 
surement. But the situation here is much more severe, since 5 g is 
known to be strongly correlated over wide range of angular scales 
and making it overwhelming the lensing signal at virtually all an- 
gular scales (Fig.[T]&[2|i. 

To a good approximation, <5 g = 6 g 5 m , where <5 m is defined 
as matter surface over-density and b g is the galaxy bias. This is 
the limiting case of deterministic bias. In general, for statistics no 
higher than second order, <5 g can be described with an extra param- 
eter, the galaxy-matter cross correlation coefficient r. b g and r are 
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Figure 1. Contaminations of the galaxy intrinsic clustering to the weak 
lensing reconstruction through cosmic magnification. The dominant con- 
tamination comes from the galaxy auto correlation, which is often sev- 
eral orders of magnitude larger than the lensing correlation. The galaxy 
intrinsic clustering also induces a galaxy-lensing cross correlation. This 
contamination can be comparable to the lensing signal. For a redshift bin 
1.0 < «b < 1-2, we plot the resulting matter power spectrum Cm b m b , 
matter- lensing cross power spectrum C mbKb and the lensing power spec- 
trum C Kh n h . For galaxies with bias b s ^ 1, a factor feg shall be applied to 



and a factor fe g shall be applied to C u 



defined through 

C s (l,z) 



bin, 



On(l,z) 



r(l, 



Cgm(l, Z) 



y/C s (l,z)C m (l,z) 



(8) 



angular power spectrum, respectively. On large scales correspond- 
ing to l^kj^x (&l ~ O.lh/Mpc), the deterministic bias (r = 1) 
is expected. However on nonlinear scales, galaxy bias is known to 
be stochastic jPenll998l;lDekel & Lahavll999l : |Pen et alj2003l : |Per] 
l2004lHoekstra et alj|2002l : lFanll2003l> . 



2.1 The estimator 

The data we have are measurements of 5i(0) at each angular pixel 
of each redshift and flux bin. Throughout the paper we use the sub- 
script "i" and "j" to denote the flux bins. For convenience, we will 
work in Fourier space. Then for a given redshift bin and a given 
multipole £, we have the Fourier transform of the galaxy over- 
density of the i-th flux bin, 5g j (£). For brevity, we simply denote it 
as hereafter. The g factor of the i-th flux bin is denoted as gi. As 
explained earlier, it is a measurable quantity, bi is the galaxy bias 
of the i-th flux bin. 

We want to find an unbiased linear estimator of the form 



k = y Wi 5. 



(9) 



Such that the expectation value of k is equal to the true k of this 
pixel. Thus the weighting function w must satisfy the following 
conditions, 



i 



o . 



(10) 



(11) 



Figure 2. Contamination of foreground galaxy intrinsic clustering to 
the cross-correlation weak lensing power spectrum between foreground 
and background redshift bins. The foreground galaxy intrinsic clus- 
tering is related with the background cosmic convergence and then 



induces a power spectrum C n 



For a couple of redshift bins 



([0.4 < Zf < 0.6, 1.0 < «b < 1-2]), we plot the weak lensing signal 
C K{ K b and compare it with this main contamination C'm { n b ■ The contam- 
ination overwhelms the signal by one or more orders of magnitude. 



Since the measured 8\ is contaminated by shot noise <5f hot , 
minimize the shot noise. This corresponds to minimize 



E 



we shall 



(12) 



Here, fn is the average galaxy surface number density of the i-th 
flux bin, a directly measurable quantity. 

The three sets of requirement (Eq. 1 1 01 [TT1 &I121 uniquely fix 
the solution. Using the Lagrangian multiplier method, we find the 
solution to be 



(Xigi + \2bi) 



(13) 



Here, the two Lagrangian multipliers Ai,2 are given by 



Ai = 
A 2 = 



2E"-Mi 

(E nihgi) 2 - E «A 2 E n>i9i ' 



(14) 



Wi is invariant under a flux independent scaling in bi . For this rea- 
son, we only need to figure out the relative flux dependence in the 
galaxy bias, instead of its absolute value. 

Despite the neat mathematical solution above, in reality we do 
not know the galaxy bias a priori. We adopt a recursive procedure 
to simultaneously solve bi, Wi and hence k. 

• The first step. In general, the power spectrum is dominated by 
the galaxy intrinsic clustering, namely 



(\5f\ 2 J = b?C mb m b + 9%giC KbKb + 2b l g l C n 



b\C„ 



^m b m b ■ 

(15) 

Here, C mb m b is the matter density angular power spectrum, C KbKb 
is the lensing power spectrum and Cm b « b is the cross-power spec- 
trum, for the background redshift bin. Hence a natural initial guess 
for bi is given by the following equation, 



(6 



(1)\2 



C n 



\sf\ 2 



(16) 



4 Yang and Zhang 



Plug b { p into in Eq.[l4]and Eq.[l3] we obtain the weighting w) ' , 
and then the first guess of lensing convergence, k' 1 ' = ^\ id^S*. 

Cm b m b is cosmology dependent, so one may think the bias re- 
construction and hence the proposed lensing reconstruction are cos- 
mology dependent. However, this is not the case. As explained ear- 
lier, the weighting function Wi does not depend on the absolute 
value of hi. So in the can be fixed to any value 

of convenience in determining the bias since it is independent with 
flux. In this sense, the bias reconstruction is cosmology indepen- 
dent and is hence free of uncertainties from cosmological parame- 
ters. 

• The second step. We subtract the lensing contribution from the 
measurement Sf, using constructed above. Our second guess 
for bi is then 



= ( \Si - go* 



(i) 



(17) 



We then obtain nxh and ft' 2 '. However, this solution is still not 



exact. The expectation value of the r.h.s is 



\8i - 9iK 



„(1) 



= ( bi - gi 



(18) 



So the expectation value of b\ is bi — gi ^ 

• The iteration. We repeat the above step to obtain and 
(p = 3, ■ • ■ ) until the iteration converges. Since we know that the 
lensing contribution is sub-dominant and we start our iteration by 
neglecting the lensing contribution, and since the flux dependence 
of the intrinsic clustering and cosmic magnification (bi and gi) are 
different , such iteration should be stable and converge. We numer- 
ically check it to be the case. 

The information of galaxy bias is not only encoded in the ob- 
served power spectrum of the same flux bin, but also encoded in 
the cross power spectra between different flux bins. In our exer- 
cise, we do not take this extra information into account. So there 
are possibilities for further improvement. 

In the above description, we have neglected shot noise (or 
equivalently assumed that it can be subtracted completely). In real- 
ity, we are only able to subtract shot noise up to the limit of cosmic 
variance. Residual shot noise introduces systematical error in the 
lensing reconstruction, which we will quantify in next section. 

Finally, we obtain the optimal estimator of cosmic conver- 
gence 



E 



(19) 



Our estimator explicitly satisfies £L w; b^ = 0. However, since 
&( can deviate from its real value bi, our estimator can be biased. 
This is yet another source of systematical errors in our weak lens- 
ing reconstruction through cosmic magnification. However, later 
we will show that such systematical error is under control. 



2.2 The reconstructed weak lensing power spectrum 

From the reconstructed k, we can reconstruct the lensing power 
spectrum. For the same redshift bin, we have 



Cbb - 

= C K 



K b I 



(20) 



E Wi 6i I C mb m b + 2 I ^ I C 



Throughout the paper we use the superscript or subscript "b" for 
quantities of background redshift bin and "f" for quantities of fore- 
ground redshift bin, and we neglect the superscript "(n)" of weight- 
ing wl from this section. For the cross correlation between fore- 
ground and background populations, we have 



Cbf 



(n) (n)< 
ft. K f 



^^w'b'j C mfKb . (21) 



Here, C Kf Kb is the lensing power spectrum between foreground and 
background redshift bins, and C mfKb is the cross power spectrum 
between foreground dark matter and background lensing conver- 



gence. We have neglected the correlation C n 



between fore- 



ground and background matter distributions. It is natural for non- 
adjacent redshift bins with separation Az > 0.1, because of phys- 
ical irrelevance. For two adjacent redshift bins, there is indeed a 
non- vanishing matter correlation. However, this correlation is also 
safely neglected since both the foreground and background intrin- 
sic clustering are sharply suppressed by factors 1/5^ V) t i ' h b\' h , re- 
spectively. 



3 STATISTICAL AND SYSTEMATICAL ERRORS 

Our reconstruction method does not rely on priors on galaxy bias, 
in this sense, it is robust. However, there is still a number of sta- 
tistical and systematical errors in the k reconstruction and weak 
lensing power spectrum reconstruction. In this section, we outline 
and quantify associated errors in the lensing power spectrum. 



3.1 The statistical error 

Our estimator minimizes the shot noise in the measurement. For 
auto power spectrum Cbb, the associated statistic measurement er- 
ror is 



(21 



v HO 2 



(2Z + 1)A;/ Sky 
For the cross-power spectrum Cbf, the statistical error is 



AC bf = 



/ '/"fshot /^shot 

(2/ + l)AZ/ sky V°b W 



(23) 



(21 + 1)AZ/ Sky ^ 



E 



K b ) 2 



E 



Throughout the paper we apply Al — 0.21, and / sky is the sky 
coverage. 

The above errors are purely statistical errors causing by shot 
noise resulting from sparse galaxy distribution. For this reason, we 
call them the weighted shot noise. We reconstruct the k from the 
measurements of lensed galaxy number over-density Sf. We solve 
the galaxy bias and weighting function from the observed galaxy 
power spectra (IfjJ'l 2 }, and this whole process is free of any given 
cosmological model. So our proposed method to reconstruct the 
weak lensing map is the one at the given sky coverage, with the 
right cosmic variance. Only when we compare our measured lens- 
ing power spectra with their ensemble average predicted by a given 
theory, we must add cosmic variance. Since statistical errors arising 
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Figure 3. The HI galaxy bias b(s, z) and the magnification bias g(s, z) = 
2(a — 1) as a function of flux s by fixing the redshift to be the central value 
of each redshift bin. The plotted curves are started from the flux limit at 
the fixed redshift. The cosmic magnification bias strongly depends on the 
flux, while the galaxy bias weakly changes with it. Such difference in flux 
dependence ensures us to find a optimal estimator to reconstruct the weak 
lensing from the cosmic magnification. 



from cosmic variance are independent to shot noise, it can be taken 
into account straightforwardly. 



3.2 Systematical errors 

We use the symbol SC to denote systematical errors in the lens- 
ing power spectrum measurement. We have identified three types 
of major systematical errors. Throughout the paper we use the su- 
perscript "(o)" (o = 1, 2, 3, • • • ) to denote them. 



tween galaxies with different flux. It is known that galaxy bias ex- 
ists a stoch astic component and hence rg < 1, especially at nonlin- 
ear sc ales dWang et al.ll2007l ; ISwanson et al„ 2008; Gil -Marin et al.l 
2010). This does not affect the determination of the galaxy bias bi, 
since we only use the auto correlation between the same flux and 
redshift bin to determine the bias. However, stochasticity does bias 
the power spectrum measurement, since now the condition (Eg. II lb 
no longer guarantees a complete removal of the galaxy intrinsic 
clustering. The systematical error induced to the auto correlation 
measurement is 



SC, 



(2) 



b b ib j b a 



(26) 



Here, Arij = 1 — r^. 

Given present poor understanding of galaxy stochasticity, we 
demonstrate this bias by adopting a very simple toy model, with 



An 





1% 



(i = i) 

(< + j) 



(27) 



This model is by no way realistic. The particular reason to choose 
this toy model is that readers can conveniently scale the resulting 
8C^ to their favorite models of galaxy stochasticity by multiply- 
ing a factor 100Ar;j (I, z). 

As we will show later, this systematical error could become 
the dominant error source. However, measuring the lensing power 
spectrum between two redshift bins can avoid this problem. Clearly, 
stochasticity in galaxy distribution does not bias such cross power 
spectrum measurement. 



3.2.1 The systematical error from deterministic bias 

The first set of systematical errors come from errors in determining 
the galaxy bias bi through Eq.[l5], even if we neglect the shot noise 
contribution to it. This bias arises due to a degeneracy that is among 
the power spectra C mb m b , C mbKb , C KbKb and galaxy bias b t in Eq. 
1151 which causes a slight deviation from its true flux dependence. 
In a companion paper we will address and clarify this issue in more 
detail (Yang & Zhang, in preparation). We will also show that this 
systematical error is correctable. 

The consequence is that Y^i w\ bi 7^ in Eq.[T9] It causes a 
systematical error in the auto correlation of the same redshift bin 

5C bb = fe^ 6 *) Cm bm b + 2 fe^l On bKb . (24) 



It also biases the cross correlation measurement between two dif- 
ferent redshift bins, 




3.2.2 The systematical error from stochastic bias 

Our reconstruction has approximated the galaxy bias to be deter- 
ministic, namely rg — 1. ry is the correlation coefficients be- 



3.2.3 The systematical error from shot noise 

By far we have neglected the influence of galaxy shot noise in de- 
termining the galaxy bias (Eq|15i. The induced error is denoted as 
8b j = b* — bj. Here, bj is the true bias and bj is the final obtained 
bias b^ from the iteration. It is reasonable to consider the case that 
shot noise is subdominant to the galaxy intrinsic clustering. Under 
this limit, 

y^shot 

(6b J )=Q,((6b j ) a ) (28) 

Om b m b 

Where Cj hot is the shot noise power spectrum of j-th flux bin in 
the observed power spectrum (Eq.ll5i. We are then able to Taylor 
expand flux weighting Wi around bj to estimate the induced bias in 
it, 

Wi(bj + Sbj) = w^ + J^^l-Sbj (29) 

j 3 3 

1 d 2 Wi I „, 
+ oZ2t^T \.Sb ] 8b k + --- 

Where Wi(bj) is the final obtained flux weighting. After a lengthy 
but straightforward derivation, we derived the induced bias in the 
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Figure 4. The contaminations before and after using the estimator to reconstruct the weak lensing auto power spectrum. We choose two redshift bins 0.4 < 
z b < 0.6 and 1.0 < z b < 1.2, for which, the upper panel shows us the suppression of matter auto power spectra, and the lower panel presents the suppression 
of matter-lensing cross power spectra, respectively. Clearly, for the same redshift bin our reconstruction method can sharply reduce the two correlations both 
induced by the intrinsic clustering. The power spectrum Cm b m b can be suppressed by an order of ~ 10 4 and the power spectrum Cm b K b can be reduced 
by one or more orders of magnitude. The suppression is more stronger for the lower redshift bin, since the value of ^ . w J? rises with redshift. Roughly 
speaking, this is caused by the increasing error in the final galaxy bias. With the increasing redshift, the contribution from the weak lensing power spectrum to 
the observed power spectrum increases, so the error in the initial galaxy bias increases and then leads to the rising error in the final galaxy bias because of the 
existing degeneracy in the process of solving galaxy bias. 



auto correlation measurement, 



3.3 Other sources of error 



sa 



(3) 



c n 



E((^ b ) 2 )(ESk) 2 



2„..b 



+(E^ b ^)(E((^ b ) 2 )E7ra & ") 



db)db) 



E«) 2 >E^|k^ 

k ij k k 



+ C n 



E((^ b ) 2 )E 



tr 



W i lb 



db)db) 1 



The induced bias in cross correlation measurement is 



7 i J J 



There are other sources of error that we will neglect in the simpli- 

(30) fied treatment presented here. First of all, we only deal with ideal- 
ized surveys with uniform survey depth without any masks. Com- 
plexities in real surveys will not only impact the estimation of errors 
listed in previous sections, but may also induce new sources of er- 
ror. These errors can be investigated with mock catalog mimicking 
real observations. We will postpone such investigation elsewhere. 

Another source of error is the determination of a (or equiv- 
alently g). For SKA that we will target at, errors in a are negli- 
gible since the galaxy luminosity function can be determined to 
high accuracy given billions of SKA galaxies with spectroscopic 
redshifts. However, this may not be the case for other surveys, due 
to at least two reasons. (1) Some surveys may not have sufficient 
galaxies at bright end. Large Poisson fluctuations then forbidden 
precision measurement of a there. (2) a is defined with respect to 

(31) the galaxy luminosity function in a given redshift bin. However, a 
large fraction of galaxy surveys may only have photometric redshift 
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Figure 5. The suppression of contamination related with foreground intrin- 
sic clustering in the reconstructed lensing cross-correlation power spectrum 
between foreground and background redshift bins. Because of the reduced 
foreground intrinsic clustering, the enormous correction C mf K b to the sig- 
nal can be suppressed by a factor ~ 10 2 . 



measurement. Errors in redshift, especially the catastrophic photo-z 
error, affect the measurement of a. 

Dust extinction is also a problem for optical surveys, but not 
for radio surveys like SKA. Dust extinction also induces fluctua- 
tions in galaxy number density, with a characteristic flux depen- 
dence a. This flux dependence differs from the a — 1 depen- 
dence in cosmic magnification and b g (s) dependence in galaxy 
bias. The minimal variance estimator can be modified such that 
WiOti = to eliminate this potential source of error, when nec- 
essary. 

In next section, we will quantify statistical and systematical 
errors listed in ^3.1l & [3T2l for the planned 21cm survey SKA. Al- 
though the estimation is done under very simplified conditions, it 
nevertheless demonstrates that these errors are likely under control. 



4 THE PERFORMANCE OF THE MINIMAL VARIANCE 
ESTIMATOR 

We target at SKA to investigate the feasibility of our proposal. SKA 
is able to detect billions of galaxies through their neutral hydrogen 
21cm emission. The survey specifications are adopted as field of 
view FOV = lOdeg 2 , total survey period tan = 5yr and total 
sky coverage 10 4 de R 2 jPewdnev et al]|2009l ; IXbdalla et alj|20ld : 
Faul kner etafl fcoiO). More details of the survey are given in the 
appendix. 

Fig.[T]and Fig. ^demonstrate contaminations of galaxy intrin- 
sic clustering to the cosmic magnification measurement from one 
same redshift bin and one couple of foreground and background 
redshift bins. For a typical redshift bin 1.0 < Zb < 1.2, the auto 
matter angular power spectrum C mb m b is larger than the lensing 
power spectrum by two or more orders of magnitude. Fig. [T]also 
shows Cm b K b is comparable to C KbKb . Since typical bias of 21cm 
galaxies is ~ 1 (Fig. [3}, this means that the galaxy intrinsic clus- 
tering overwhelms the lensing signal by orders of magnitude. Sim- 
ilarly, in Fig. [2] the cross power spectrum C mfKb induced by the 
foreground intrinsic clustering overwhelms the weak lensing power 
spectrum C KfKb by one or more orders of magnitude. These big 
contaminations related galaxy intrinsic clustering make the weak 
lensing measurement difficult from the directly cosmic magnifi- 
cation measurement, unless for sufficiently bright foreground and 



backg round galaxies at sufficiently high redshifts (Zhang "& Penl 
l200a> . 

As we explained in earlier sections, the key to extract the 
lensing signal from the overwhelming noise is the different depen- 
dences of signal and noise on the galaxy flux. Fig. [3] shows the 
lensing signal and the intrinsic clustering indeed have very differ- 
ent dependences on the galaxy flux. For the 21cm emitting galax- 
ies, the flux dependence in g = 2(q — 1) is much stronger than 
that in the galaxy bias. Furthermore, g changes sign from faint end 
to bright end. Such behavior can not be mimicked by bias, which 
keeps positive. 

From such difference in the flux dependences, we expect 
that our estimator to significantly reduce contaminations from the 
galaxy intrinsic clustering. As explained earlier, in the galaxy cor- 
relation between the same redshift bin, the intrinsic clustering in- 
duces a systematical error proportional to C mb m b (Eq.l24t and an 
error proportional to C mbKb (Eq.l24t. In the ideal case that both the 
stochasticity and shot noise in galaxy distribution can be neglected, 
the systematical error proportional to C mb m b will be suppressed 
by a factor ^b^) 2 (Eq.|24j. Fig.|4]shows that this suppres- 

sion factor is of the order ~ 10 4 at interesting scales 10 < £ < 10 4 . 
For the same reason, the systematical error proportional to C mb K b 
will be suppressed by a factor 1/ £\ w°b^ ~ 10 2 over the same 
angular scales. The systematical error induced by foreground in- 
trinsic clustering in the weak lensing reconstruction between two 
redshift bins is proportional to C mf Kb (Eq.|25ll. Our estimator also 
suppresses it by a factor 1/ X^i w \b\ ~ 10 2 , as shown in Fig.[5] 



4.1 Same redshift bins 

The lensing power spectrum can then be directly measured through 
the reconstructed lensing maps. This can be done for maps of the 
same redshift bin. Fig.[6]compares the residual systematical errors 
to the lensing power spectrum signal. Overall, our minimal vari- 
ance estimator significantly suppresses the galaxy intrinsic cluster- 
ing and makes the weak lensing reconstruction feasible. For source 
redshift Zb ~ 1 (e.g. 1.0 < Zb < 1.4), all investigated systematical 
errors and statistical errors are suppressed to be subdominant to the 
signal. We expect the lensing power spectrum to be measured with 
an accuracy of ~ 10%-20%. 

However, at lower and higher redshifts, the reconstruction 
is not successful. We observe the systematical error fiC^^ in- 
creases with increasing redshift and overwhelms the lensing signal 
at Zb ^ 2. In the current formalism it is difficult to explain this be- 
havior straightforwardly. But roughly speaking, this is caused by 
worse initial guess on the flux dependence of galaxy bias coupled 
with the degeneracy explained earlier. The initial guess is accurate 
to the level of C KbKb /C mb m b . So the associated error increases 
with redshift. 

This systematical error looks rather frustrating. However in 
the companion paper (Yang & Zhang, in preparation) we will 
show that the systematical error 8C^ is correctable. Where 
we can separate the degeneracy, and reconstruct a quantity 
2/» b = \/Cm b m b (b* + giC mbKb /C mb m b ) thi'ough a direct multi- 
parameter fitting against the measured power spectra, which per- 
fectly mimics the flux dependence of galaxy bias because the 
power spectrum C mb m b is independent with the flux and the cor- 
rection term C mbKb /C* mb m b is small especially for high redshift. 
Although it is bad to find that the final convergence depends on the 
initial guess of galaxy bias, y° as a guess of galaxy bias does re- 
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Figure 6. Composition of the reconstructed weak lensing power spectrum for the same redshift bin. We plot the weak lensing power spectrum C KbKb by bold 
solid line. The three types of systematical errors 

SC bb from 

error in deterministic galaxy bias, SC^ from the stochastic bias, and SC^ from the shot noise 
are presented by the solid, dotted and dash-dot-dot-dotted lines, respectively. While statistical error is plotted by the dashed line. Here the statistical error is 
called weighted shot noise only from the sparse galaxy distribution, since we aim to reconstruct the weak lensing at each angular pixel with the corresponding 
cosmic variance. For the intermediate redshift bins 1.0 < Zb < 1-2 and 1.2 < Zb < 1-4, the signal overwhelms all these errors which can be controlled to 
~ 10%-20% level. 



duce the systematical error 5C' 1 ' and then works far better than the 
obtained galaxy bias in this paper. 

The systematical error arising from shot noise becomes 

non-negligible at z\, ~ 2, due to sharply decreasing galaxy den- 
sity and increasing shot noise at these redshifts (see Fig. lA2h . We 
find that this error is always subdominant to either 5C bt ' or 5C{^ . 
Interestingly, both SC^j and the weighted shot noise ACbb (sta- 
tistical error, Eq. 1221 have similar shapes at all redshifts. Further- 
more, both roughly scale as 1° at low redshifts. These are not coin- 
cidences. (1) The term oc C mbmb is dominant in SC^ . This term 
is also oc ((8b];) 2 } oc l/C mbmb /fij. So both 5C bb > and AC bb 
are the sums of weighted in different ways and hence have 
similar shapes. (2) The galaxy bias in our fiducial model is scale 
independent. This results in scale independent weighting function 
Wj, as long the bias can be determined to high accuracy. For these 
reasons, SC^j , ACbb oc 1°. This is the case at low redshift. (3) 
However, the accuracy in determining galaxy bias is significantly 
degraded by contamination proportional to C KbKb /C mb ,n b , which 



is scale dependent and increases with redshift. This is the reason 
both 5Cbt, and ACbb show complicated angular dependence at 
z h ~ 2 (Fig.©. 



At low redshift (e.g. 0.4 < Zb < 0.6), the systematical er- 
ror <5Ckb' arising from galaxy bias stochasticity becomes domi- 
nant or even exceeds the lensing signal. This is what expected. 
The galaxy intrinsic clustering is stronger at lower redshift. This 
amplifies the impact of galaxy stochasticity. Weaker lensing signal 
at lower redshift further amplifies its relative impact. As expected, 
Fig.|6]shows that SC^j /C KbKb decreases with increasing redshift 
and becomes negligible at Zb ~ 2. Since <5C bb ' oc Arij, the impor- 
tance of <5C bb ' is sensitive to the true nature of galaxy stochastic- 
ity. Its value should be multiplied by a factor lOOArij for fiducial 
value of Ary 7^ 1%. We hence conclude that galaxy stochasticity 
is likely the dominant source of error in weak lensing reconstruc- 
tion through cosmic magnification. 
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Figure 7. The weak lensing signal, the systematical error and the statistical error in the weak lensing reconstruction from the foreground and background 
redshift bins. The solid line is the cross-correlation power spectrum of cosmic convergence. The dotted line represents the systematical error combing three 



types <5Cbf = fiCyrf + 5C¥f + <5C b f ; and the dashed line corresponds to the weighted shot noise. In <5Cbf, SC^ 1 from the deterministic galaxy bias is 



-.(-') 



-.(3) 



(1) 



dominant and the systematical error from the stochastic bias can be avoided in such cross lensing power spectrum, namely <5C b j = 0. For every couple of 
foreground and background redshift bins, the cross reconstructed weak lensing signal dominates at scale range 10 < I < 10 4 and it can be measured to reach 
~ 10%-20% accuracy. 



4.2 Different redshift bins 

Fortunately this stochasticity issue can be safely overcome in the 
lensing power spectrum measurement through lensing maps recon- 
structed in two different redshift bins (foreground and background 
bins). The results are shown in Fig. [7] In this case, the systemati- 
cal error is dominated by SCPf . The stochasticity does not induce 

(2) 

systematical error so that 6C£' = 0. 

Overall, the lensing power spectrum measurement through 
cross correlating reconstructed maps of different redshift bins is 
more robust than the one based on the same redshift bin. The re- 
construction accuracy can be controlled to 10%-20% over a wide 
range of foreground and background redshifts. 

At last, for a consistency test, we check whether our results 
depend on the division of the flux bin. As expected, various sys- 
tematical errors and statistical error change little with respect to 
different choices of flux bins, as long as these bins are sufficiently 
fine. 



4.3 Uncertainties in the forecast 

There are a number of uncertainties in the forecast, besides the ones 
discussed in previous sections. (1) In the fiducial galaxy intrinsic 
clustering model, we have neglected the scale dependence in galaxy 
bias. (2) We have neglected cosmic variance in the lensing signal 
so the fiducial power spectrum is the ensemble average. But we do 
not expect it can significantly impact our result, since the cosmic 
variance at most relevant scale is small due to ZAZ/ S k y 2> 1 with 
a given large sky coverage of SKA. (3) The toy model of galaxy 
stochasticity is too simplified. In reality, the cross correlation coef- 
ficient r should be a function of redshift, angular scale and galaxy 
type and flux. 

For these reasons, the numerical results presented in this pa- 
per should only be trusted as rough estimation. Robust evaluation of 
the weak lensing reconstruction performance through cosmic mag- 
nification requires much more comprehensive investigation. Nev- 
ertheless, the concept study shown in this paper demonstrate that 



10 Yang and Zhang 



weak lensing reconstruction through cosmic magnification is in- 
deed promising. 



5 CONCLUSIONS AND DISCUSSIONS 

We propose a minimal variance estimator to reconstruct the weak 
lensing convergence k field through the cosmic magnification ef- 
fect in the observed galaxy number density distribution. This es- 
timator separates the galaxy intrinsic clustering from the lensing 
signal due to their distinctive dependences on the galaxy flux. Us- 
ing SKA as an example, we demonstrate the applicability of our 
method, under highly simplified conditions. It is indeed able to sig- 
nificantly reduce systematical errors. We have identified and quan- 
tified residual systematical errors and found them in general under 
control. Extensive efforts shall be made to test our reconstruction 
method more robustly and to improve this method. 

Comparing to previous works, our method has several key fea- 
ture s/ advantages . 

• Unlike weak lensing reconstruction through cosmic shear, it 
does not involve galaxy shape measurement and reconstruction and 
hence avoids all potential problems associated with galaxy shapes. 
Hence the reconstructed lensing maps provide useful independent 
check against cosmic shear measurement. 

• Unlike existing cosmic magnification measurements which 
actually measure the galaxy-lensing cross correlation, our estima- 
tor allows directly reconstruction of the weak lensing k field. From 
the reconstructed k, we are able to directly measure the lensing 
power spectra of the same source redshift bin and between two red- 
shift bins. These statistics do not involve galaxy bias, making them 
more robust cosmological probes. The usual lensing tomography is 
also directly applicable. 

• Unlike our previous works dZhang & Pen] I2005L 120061) . the 
new method does not require priors on the galaxy bias, especially 
its flux dependence. Our methods is able to simultaneously mea- 
sure the galaxy bias (scaled with a flux independent factor) and the 
lensing signal. Hence we do not adopt any priors on the galaxy 
bias (other than that it is deterministic) and treat the galaxy bias as 
a free function of scale and flux. The price to pay is degradation 
in constraining galaxy bias and in lensing reconstruction. Adding 
priors on the galaxy bias can further improve the reconstruction 
precision, although the reconstruction accuracy will be affected by 
uncertainties/biases in the galaxy bias prior. For this reason, we do 
not attempt to add priors on galaxy bias in the reconstruction. 

• Our method is comp lementary to a recent proposal by 
iHeavens & Joachimil d201ll ). which proposes a nulling technique 
to reduce the galaxy intrinsic clustering by proper weighting in 
redshift. Comparing to this method, our method only utilizes extra 
information encoded in the flux dependence to reduce/remove the 
galaxy intrinsic clustering. It keeps the cosmological information 
encoded in the redshift dependence disentangled from the process 
of removing the intrinsic clustering. 

The proposed approach is not the only way for weak lens- 
ing reconstruction through cosmic magnification. The current pa- 
per focuses on direct reconstruction of the lensing convergence k 
map. In a companion paper, we will focus on direct reconstruc- 
tion of the lensing power spectrum (Yang & Zhang, in prepara- 
tion). We will show that combining two-point correlation measure- 
ments between all flux bins, the lensing power spectrum can be 
reconstructed free of assumptions on the galaxy intrinsic cluster- 
ing. We will see that this approach is more straightforward, more 



consistent and easier to carry out. However the method presented 
in this paper does have advantages. Since it reconstructs the lensing 
k map, higher order lensing statistics such as the bispectrum can be 
measured straightforwardly. Furthermore, the reconstructed k map 
can be straightforwardly correlated with other tracers of the large 
scale structure. For example, it can be correlated with the lensing 
map reconstruct ed from CMB lensing dSeliak & Zaldarriagall 199911 
Hu & Okamotdl2002h or 21cm lensing. Furthermore, through this 
approach we can have better understanding on the origin of vari- 
ous systematical errors, which can be entangled in the alternative 
approach. 

Our reconstruction method is versatile to include other com- 
ponents of fluctuation in the galaxy number density. The extinction 
induced fluctuation discussed earlier is one. High order corrections 
to the cosmic magnification is another. Taylor expanding Eq. Q]to 
the second order, we obtain 



5g = 8 3 + 2(q - 1)k 
+2(q-1) 



(32) 



(k5 b - (nSg)) + ~( 7 



<7 2 » 



+ (1 -5a + 2g 2 ) [k 2 - ( K 2 )] + O [k\ 7 3 , ■ • ■ ] . 

Here g 2 = (s 2 /n)d 2 n/ds 2 is related with the second derivative of 
luminosity function. 

The above result shows that the k reconstruction through cos- 
mic magnification is biased by terms proportional to n8 g and j 2 
(second line in the above equation). Similar biases also exist in 
cosmic shear measurement. We recognize nS g as the source-lens 
coupling. The -y 2 term is analogous to the K-y term caused by re- 
duced shear 7/(1 — ft). Precision lensing cosmology has to model 
these corrections appropriately. 

The high order corrections oc 1 — 5a + 2g 2 can in principle be 
separated due to its unique flux dependence. However, it is unclear 
whether the reconstruction is doable, even for a survey as advanced 
as SKA. 
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APPENDIX A: SKA SURVEY 

SKA is a future radio survey with aiming to construct the world's 
largest radio telescope. Through the neutral hydrogen emitting 21 
cm hyperfine transition line, it can observe large sample of 21 cm 
galaxies. Even at high-redshift, the observed 21 cm galaxies is ex- 
tremely excess than the QSOs or LBGs that is usually used as back- 
ground object dScranton et alll2005t lHildebrandt ~etai]|2009h , so it 
is easily to overcome the shot noise in the measurement of cosmic 
magnification effect, of which precise measurement is sensitive to 
the shot noise. In addition, compared with photometric survey, the 
spectroscopic survey can determine the redshift more precise by 
using the redshifted wavelength of 21 cm line (A = Ao(l + z)). 
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Figure Al. The redshift evolutions of HI mass limit (dashed line) and the 
characteristic mass in model C (dotted line). The characteristic mass in non- 
evolution model is plotted by the solid line. 



Furthermore, radio survey is free of galactic dust extinction, which 
is correlated with foreground galaxies and then induces a correc- 
tion to the galaxy-galaxy cross correlation. As a consequence, it is 
expected a good performance of our proposed method to measure 
the cosmic magnification for SKA. In this paper, we adopt the sur- 
vey specifications as follows: the telescope field-of-view (FOV) is 
lOdeg 2 without evolution, the overall survey time is 5yr and the 



sky coverage is lOOOOdeg 



Al HI mass limit 



By assuming a flat profile for the emission line with frequency, we 
obtain the following equation linking Mm with the observed flu x 
density s by atomic Physics (detail in Abdalla & Rawlingdj2005i) ). 



Mm = 



16tt 



3 Ai\hc 



X 2 (z)sV(z)(l + z) 



(Al) 



Where x( z ) ls the comoving angular distance, A21 is spontaneous 
transition rate, tuh is the atomic mass of hydrogen, h is Planck's 
constant and V(z) is the line-of-sight velocity spread. We assume a 
scaling with redshift for the typical rest-frame velocity width of the 
line V(z) = Vb/vTTz, where V = 300kms~\ and ignore the 
effects of inclination. The r.m.s. sensitivity for a dual-polarization 
radio receiver at system temperature T sys for an integration of du- 
ration t on a telescope of effective collecting area A e g is given by 



V2T si 



k r 



VcA cf f 



(A2) 



Here, the correlation efficiency r\ c is adopted as r\ c = 0.9, the du- 
ration time t = 40/i for each area of the sky and the width of the 
HI emission line determines the relevant frequency bandwidth Ais, 
which is related to a line-of-sight velocity width V(z) at redshift z 



Av 



vo V{z) 

1 + 2 C 



(A3) 



The flux detection limit sn m for galaxy is defined by the threshold 
parameter n CT = sn m /s rms . Here, we apply n a = 5. From Eq . 
dA H . we can obtain the HI mass limit. 
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A2 HI mass function 

we assume that the HI mass function at all redshifts is described by 
a Schechter function, 



<p(Mm,z)dM m = <t> 



. /MhiV 



exp 



A/hi \ dMm 
~ M* ) M* 



(A4) 



Here, the parameter /3 is low-mass slope, M* is characteristic 
mass and <f>* is normalization. The HIPASS survey reported the 
results: /3 = -1.3, = 0) = 3.47h~ 2 lO 9 M and </>* = 

0.0204h 3 Mpc" 3 (Zwaan et al. 2003). There is little solid measure- 
ment of <f>* (z) and M* (2) other than at local universe. But for this 
form of mass function, there exists a tight relation between Sim, 
the present day critical density p c (z = 0) and <f>* (z)M* (z), 



Slmh = r(/3 + 2)0* M* (z) I 'p c (z = 0) 



(A5) 



Where V is the Gamma function. Observation of damped Lyman- 
alpha (DLA) system and Lyman-alpha l imit s ystem can be used to 
measure Sim dPeroux et al.ll200lLl2003ll2004h . We use the follow- 
ing fu nctional form produc ed by fitting to DLA data points used in 
paper dAbdalla 



orm produo 
et al.ll2010h 



Sli 



N [1.813 — 1.473(1 + z) 



(A6) 



The normalization constant N is fixed by the value of <j)*(z = 

0)M*(* = 0). 

Here, we adop t model C in papers dAbdalla & Rawlings 2005 ; 
Abdal la et al j201Ch as evolution model of HI mass function. In this 
model, cj>* scales with z by using the DLA results. The break in the 
mass function M* is controlled by the cosmic star formation rate 
and given by 



M*(z) 



Sl a t a ,r(0)/Slm(0) + 2\ d3 



M*(0) \Sl atar (z)/Sl m {z) + 2 



(A7) 



Where it is assumed that Slm(z) = Sln 2 ( z ) and D(z) is the growth 
factor. The fractional density in stars Q s ta.i(z) = Pstm(z)/ p c (z = 
0) is deduced by using the fits to the c osmic star-formation history 
in lChoudhurv & P admanabhan] d2002h . 



Sl a 



SFR(>') 



p e (* = 0)J a H(z')(l + z') 



-dz 



and 



SFR(z) = 



0.13 



1 + 6exp(-2.2z) 



0.5 



1 



Xtid(z) 



(A8) 



(A9) 



Here, x( z ) an d Xfid{z) are the co-moving distances at redshift 2 
for the adopt cosmology model and a fiducial cosmology model 
with SIm — 1 and SIa ~ 0, respectively. The units of SFR 
are h 2 M0yr -1 Mpc -3 . Fig . I All shows us the evolutions of HI 
mass limit and characteristic mass in evolution model C. The no- 
evolution characteristic mass is also plotted for comparison. In Fig . 
IA2I we plot the redshift distributions of HI galaxy for the evolution 
model C and the no-evolution model. 




Figure A2. The redshift distributions of HI galaxies for no-evolution model 
and evolution model C. 



A4 The galaxy bias b a for HI galaxy 

In our forecast, we adopt a deterministic bias b by assuming the 
cross correlation coefficient r — 1, and hence we estimate the error 
induced by this assumption. 

The large s cale bias o f HI galaxies estimated by method pre- 
dicted in paper 



b(M r 



,0.5 , 



(0.06-0.02n off ) 



1 + 



(All) 



Where v — 5 c (z) /<t(Mdm). c(Mdm) is the linearly evolved rms 
density fluctuation with a top-hat window function and S c (z) = 
1.686/D(2). The effective index is defined as the slope of P(k) at 
the halo mass Mum, 
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HI galaxies are selected by their neutral hydrogen mass Mm- To 
calculate the associated biases of these galaxies, we need to convert 
Mm to the galaxy total mass Mum- Since most baryons and dark 
matter are not in galaxies while most neutral hydrogen atoms are in 
galaxies, we expect that /hi 3> Slm/Sl m . Here, we choose /hi = 
0.1. Although there are some problems to model the HI galaxy bias 
by this way, it is emphasized that our reconstruction method is not 
limited to the specific value of galaxy bias adopted. As long as the 
galaxy bias has different flux dependence with the magnification 
bias. Their different flux dependences are presented in Fig . [3] 



A3 The slope of the HI galaxy number count 

In Eq .[5] n is the mass function of HI galaxies. We can derive the 
parameter a as a function of flux s and redshift z, which is 

a (Mm,z) = -/3-l+ 1i ^fir =a(s,z). (A10) 
M (z) 

The relationship between flux s and Mhi is given by Eq . IA1I 



